接触 热 阻 对 燃料 电池 温度 分 布 影响 的 数值 模拟 


曹 涛 锋 ， 母 玉 同 ， 丁 靖 ， 何 雅 玲 ， 陶 文 狂 
(西安 交通 大 学 能 源 与 动力 工程 学 院 ， 热 流 科学 与 工程 教育 部 重点 实验 室 ， 西 安 ，710049) 

摘要 : 本 文采 用 三 维 、 非 等 温 、 以 及 考虑 各 向 异性 扩散 层 性 质 的 燃料 电池 模型 分 析 了 气体 扩散 层 和 流 场 板 肋 板 之 间 的 
接触 热 阻 对 电池 温度 分 布 和 性 能 的 影响 。 计 算 结 果 表 明 ， 接 触 热 阻 对 电池 最 高 温度 及 温度 分 布 均 有 很 大 影响 。 电 池 输 
出 电压 为 0.3V 时 ， 有 接触 热 阻 工 况 电 池 最 大 温度 高 于 无 接触 热 阻 工 况 约 8.5K。 此 外 ， 有 接触 热 阻 计算 工 况 MEA 内 温 
度 分 布 近似 为 “1%” 型 ， 而 无 接触 热 阻 计 算 工 况 MEA 内 温度 分 布 近似 为 “A” 型 。 因 此 ， 对 于 燃料 电池 温度 分 布 的 准 
确 预测 ， 不 能 忽略 接触 热 阻 的 影响 。 
关键 词 : 质子 交换 膜 燃 料 电池 ， 数 值 模拟 ， 温 度 分 布 ， 接 触 热 阻 
中 图 分 类 号 : TK121 文献 标识 码 A 


Numer ical lnvestigation the Temperature Distribution of PEM Fuel Cell With 


the Effect of Thermal Contact Resistance 


CAO Tao-Feng, MU Yu-Tong, DING Jing, HE Ya-Ling, TAO Wen-Quan 
(School of Energy and Power Engineering, Key Laboratory of Thermo-Fluid Science and Engineering of MOE, Xi’an Jiao tong 
University, Xi’an, 710049, China) 
Abstract: In this paper, a three dimensional non-isothermal PEM fuel cell model with anisotropic gas 
diffusion layer(GDL) is applied to investigate the effect of the thermal contact resistance(TCR) between 
flow plate rib and GDL on cell temperature distribution and performance. From the numerical results, it is 
found that both the maximum cell temperature and the temperature distribution are greatly affected by TCR. 
When the output cell voltage equals to 0.3V, the maximum cell temperature with TCR is about 8.5K larger 
than that of without TCR. Furthermore, the temperature profile of thermal contact resistance case looks like 
a "\" style, however, the temperature distribution without thermal contact resistance looks like a "A" style. 
Therefore, it is concluded that, to accurately predict the temperature distribution of a PEM fuel cell, the 
thermal contact resistance between flow plate rib and gas diffusion layer cannot be ignored. 
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过 去 十 年 里 ， 学 者 们 对 电池 温度 分 布 进行 了 ) 
泛 研 究 。VieIl，Wang 叫 ，LinBl]， 以 及 Thoms 中 等 
操作 温度 是 影响 质子 交换 膜 燃料 电池 性 能 的 关 ”采用 实验 方法 研究 了 不 同 工 况 下 燃料 电池 的 温度 分 
键 参数 之 一 ， 为 了 确保 燃料 电池 的 稳定 运行 ， 必 须 布 ， 与 此 同时 ，JuE 4，Suil 9，Wug 等 人 则 以 数值 
将 其 工作 温度 控制 在 合适 的 范围 之 内 。 温 度 过 高 会 。” 方法 模拟 了 燃料 电池 内 的 温度 分 布 。 与 实验 方法 相 
加 快 质 子 交 换 膜 及 众 化 剂 的 衰减 ， 影 响 电 池 寿 命 ; 比 ， 数 值 模拟 方法 具有 经 济 实惠 ， 方 便 易 用 ， 可 以 
温度 过 低 会 使 电化 学 反应 速率 下 降 ， 或 者 引起 电极 得 到 详细 的 场 分 布 等 优点 。 

“水 渡 ” 造成 电池 性 能 下 降 。 因 此 ， 合 理 的 电池 热 采用 数值 模拟 方法 预测 燃料 电池 温度 分 布 时 ， 
管理 方法 对 燃料 电池 的 稳定 运行 极为 重要 。 而 理解 其 准确 程度 依赖 于 以 下 几 个 问题 : (1) 电池 热 边 界 
燃料 电池 内 传 热机 理 以 及 电池 温度 分 布 则 是 制定 电 条 件 ; (2) 电池 内 热源 ; (3) 电池 各 组 件 的 导热 系 
池 热 管理 策略 的 基本 依据 。 数 ; 〈4) 电池 不 同 组 件 之 间 的 接触 热 阻 。 然 而 ， 绝 
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大 多 数 数值 计算 模型 忽略 了 接触 热 阻 的 影响 。2005 
年 ，Birgerssont "等 人 采用 二 维 模 型 分 析 了 接触 热 阻 
对 电池 温度 分 布 的 影响 。 近 来 ，Sadeghin' 2 等 人 发 
现 接触 热 阻 是 电池 总 热 阻 的 主要 部 分 ， 忽 略 接触 热 
阻 将 会 导致 电池 温度 分 布 预测 错误 。 
基于 以 上 分 析 ， 本 文采 用 燃料 电池 两 相 非 等 温 
模型 研究 了 电池 流 场 板 肋 板 与 气体 扩散 层 间 的 接触 
热 阻 对 电池 温度 分 布 的 影响 。 以 下 分 别 介绍 了 该 燃 
料 电 池 两 相 非 等 温 模型 及 接触 热 阻 的 实施 方法 ， 之 
后 分 析 讨 论 了 接触 热 阻 对 电池 温度 分 布 及 性 能 的 影 
响 ， 最 后 给 出 了 结论 。 
1 计算 模型 与 数值 方法 
1.1 燃料 电池 两 相 非 等 温 模 型 

本 次 计算 采用 的 燃料 电池 两 相 非 等 温 模型 以 本 
课题 组 开发 的 燃料 电池 宏观 尺度 模型 为 基础 羽 ， 
加 入 了 各 向 异性 气体 扩散 层 对 燃料 电池 性 能 的 影 
响 。 模 型 假设 电池 稳 态 运行 ， 反 应 气体 为 理想 气体 ， 
流动 状态 为 层 流 ， 气 体 通道 中 无 液态 水 ， 并 忽略 电 
池 不 同 组 件 间 的 接触 电阻 。 模 型 控制 方程 包括 质量 、 
动量 、 能 量 、 组 分 以 及 电荷 守恒 方程 ， 各 控制 方程 
表达 式 见 表 1， 控 制 方程 源 项 表达 式 见 表 2， 模 型 参 
数 及 其 他 附加 方程 详 见 文献 [1$-18]。 
1. 2 计算 区 域 与 边界 条 件 
1.2.1 计算 区 域 

本 文 计算 区 域 如 图 1 所 示 ， 为 平行 流 场 板 典 型 
单元 。 考 虑 到 平行 流 场 板 结构 沿 y 方 向 对 称 ， 因 此 ， 
计算 区 域 取 为 通道 和 肋 板 的 一 半 ， 包 括 燃 料 电池 的 
所 有 构成 部 分 。 图 1 中 ， 截 面 1 为 流动 方向 中 心 截 
面 ，Line-1 为 截面 1 上 肋 板 与 通道 交界 线 。 
1. 2. 2 边界 条 件 

不 同 界面 的 边界 条 件 可 表述 为 ; 
(1) 进口 边界 : 

给 定 速 度 及 反应 气体 质量 分 数 及 温度 ， 


Urin 一 C 4 (1) 
k 3 


“3 
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(2) 出 口 边界 : 
假设 出 口 边 界 为 充分 发 展 条 件 ， 即 
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el (3) 
(3) 左右 边界 : 
假设 左右 边界 为 对 称 边界 ， 即 ， 
撞 w 2? 厂 
v= 0, 0, 0， 0 4 
挡 y ?7 ® 
(4) 上 下 边界 : 


计算 区 域 上 下 边界 与 外 部 环境 接触 ， 给 定 电池 
操作 电压 以 及 温度 ， 即 ， 


办。 本 0, 人 办。 至 Vo, T a To (5) 


计算 区 域 几何 参数 及 边界 条 件 相 关 取 值 见 表 3。 


aa 氮气 
一 > 氧气 


1 计算 区 域 
Fig.1 Computational domain 
1.3 接触 热 阻 的 实施 方法 
气体 扩散 层 一 般 由 碳 纸 制 成 ， 其 表面 粗糙 ， 无 


法 与 气体 通道 肋 板 充分 接触 ， 通 常 在 二 者 接触 面 上 
产生 接触 热 阻 ， 增 加 传 热 阻力 。 由 于 接触 界面 的 间 
孙 很 小 ， 无 法 直接 建 模 。 因 此 ， 本 文 在 实施 接触 热 
阻 时 ， 将 流 场 板 肋 板 划分 为 两 部 分 ， 如 图 1 所 示 。 
通过 改变 靠近 气体 扩散 层 部 分 肋 板 的 导热 系数 ， 使 
其 产生 的 热 阻 与 接触 面 形成 的 接触 热 蛆 相当 ， 以 达 
到 接触 热 阻 的 效果 。 单 位 面积 热 阻 通常 表示 为 ; 
d 

R= (0) 
因此 ， 靠 近 气体 扩散 层 部 分 肋 板 薄 层 的 导热 系数 为， 


di; 
1 三 Rib_S (7) 
Rj 


表 1 控制 方程 


Table 1 Governing equations 


控制 方程 


质量 守恒 方程 vtLpsu)=sS， 
动量 守恒 方程 vr pee -wp, + VI eVus) +S, 
(E2 ) £ 
组 分 守恒 方程 Vpens0) = Vip Dar VO) +S, 
电荷 守恒 方程 Vi(oVG)=5;, V:(0, VG )=5, 
能 量 守 恒 方程 VI pcpu,T) =VLKarY7)+Sr 
液态 水 饱和 度 方程 中 pn =VI(PD.Vs)-—S, 
膜 水 含量 方程 Y(D.ycj-Y mE] -5 
表 2 不同 区 域 源 项 表达 式 
Table 2 Source terms in different regions 
控制 方程 源 项 
质量 守恒 方程 崔 化 层 : 5S =S +S +S， 
动量 守恒 方程 扩散 层 及 催化 层 : S, 三 一 [aus / KK,, 
组 分 守恒 方程 阴极 催化 层 : 5,=-(Gi./4F)M,， 阳 极 众 化 层 : 5, = -Gi /2F)M, 
> 催化 层 ，S. = 5 + 3 ， 扩散 层 ， 5S, -8 
共 守 恒 方 程 阳极 催化 层 : 8% .= 一 ， 阴 极 催化 层 ，S = 
人 阳极 催化 层 ，S， =i， 阴 极 催化 层 ， Ss,, = 
能 量 守恒 方程 催化 层 :5， =i0+T E+ 人 ， 质 子 交换 膜 5， -二 
液态 水 饱和 度 方 程 S$) = do me (1—s)(p, — pi )q + A Te s(p,— pa)(—9) 
阳极 催化 层 : S =- 思 (cv -cv 
0 阴极 催化 层 ，S = -yc -csi2 
表 3 模型 几何 尺寸 及 操作 参数 流 道 下 温度 均 高 于 肋 板 下 温度 ， 这 是 因为 通道 中 气 
Table3 Model geometric dimension and operating parameters 体 导 热 系 数 远 小 于 肋 板 导热 系数 ， 肋 板 下 的 热量 较 
参数 符号 ”单位 值 通道 下 更 容易 散 出 。 有 接触 热 阻 时 ， 电 池 最 高 温度 
二 0 高 出 无 接触 热 阻 工 况 约 3C， 且 电池 中 心 高 温 区 较 
en 5p。 Pa 101325101325 ”大 ， 这 是 因为 接触 热 阻 增 大 了 传 热 阻力 ， 使 得 电池 
阳极 /阴极 化 学 计量 比 /6 33 内 产生 的 热量 个 易 敬 出 。 _ 
入 口 氧气 /空气 相对 湿度 RH RH 60%/60% 不 同 电压 下 ， 沿 流动 方向 中 心 截 面 位 置 ， 流 道 
ae a 以 看 出 ， 随 着 电压 的 升 高 ， 电 池 最 高 温度 降低 ， 这 
云 厅 扩 把 ， ss VA = 
膜 厚度 记 本 5 是 因为 ， 电 压 升 高 时 ， 电 流 减 小 电池 产生 的 热量 减 
催化 层 厚度 刀 es 0.0284 小 。 | 比 图 4(a) 和 (b) 可 以 发 现 ， 有 接触 热 阻 时 电池 
接触 热 阻 Rs mkw! 2.5x104 最 高 温度 高 于 无 接触 热 阻 时 ， 电 池 电 压 为 0.3V 时 ， 
有 接触 热 阻 时 电池 最 高 温度 约 为 365K， 而 无 接触 热 
本 阻 时 ， 电 池 最 高 温度 仅 为 356.5K， 二 者 相差 8.5K。 
2 生来 和 时 论 这 同样 是 由 于 接触 执 阻 增 大 了 传 热 胃 力 ， 阻 碍 电池 
计算 用 到 的 电池 操作 条 件 见 表 3。 为 了 分 析 接 触 妾 让 的 执 量 向 从 其 兴 二 说 下 泪 唐 介 
内 产生 的 热量 向 外 散发 。 此 外 ， 两 种 工 况 下 温度 分 


热 阻 对 电池 温度 分 布 和 性 能 的 影响 ， 本 节 对 比 了 有 
〗 布 曲 线 亦 有 很 大 不 同 。 有 接触 热 阻 时 ， 肋 板 上 温度 
接触 热 阻 和 无 接触 热 阻 两 种 工 况 的 温度 分 布 ， 阴 极 。。 几乎 不 变 ， 到 达 气 体 扩散 屋 时 ， 突 然 四 计 


闹 太 水 猩 和 诺 、 » 
ee MEA 内 温度 变化 相对 平缓 ， 尤 其 气体 扩散 层 内 温度 
.1 温度 变化 较 小 ， 温 度 分 布 近似 呈 "型 。 


电池 输出 电压 为 0.6V 时 ， 两 种 工 况 下 的 电池 温 
度 分 布 如 图 3 所 示 。 可 以 发 现 ， 两 种 工 况 下 ， 电 池 


| 


3 温度 分 布 (a) 接触 热 阻 工 况 (b) 无 接触 热 阻 工 况 
Fig.3 Temperature distribution (a) with TCR (b) without 
TCR 
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图 4 不 同 电压 下 Line-1 位 置 温度 分 布 曲线 (a) 接触 热 
阻 工 况 ，(b) 无 接触 热 阻 工 况 

Fig.4 Temperature profiles along Line-1 under different cell 


voltage (a) with TCR, (b) without TCR 
而 无 接触 热 阻 时 ， 肋 板 内 温度 呈 线 性 变化 ， 从 
ed 而 MEA 内 温度 变化 较 快 ， 以 阴极 
炭化 层 为 中 心 ， 近 似 呈 “A” 型 分 布 。 因 此 ， 采 用 数 
ee a 
昌 的 影响 。 
2.2 阴极 液态 水 饱和 度 分 布 
电池 电压 为 0.6V 时 ， 流 动 方向 中 心 截 


i 


i 


TS 


处 ， 阴 
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极 液态 水 无 量 纲 饱和 度 分 布 如 图 5 所 示 。 由 于 饱和 
蒸汽 压 受 慢 度 影响 较 大 ， 温 度 高 时 ， 水 的 饱和 蒸汽 
压 升 高 ， 饱 和 度 降低 ， 温 度 低 时 ， 水 的 饱和 蒸汽 压 
降低 ， 饱 和 度 升 高 。 因 此 ， 电 池 肋 板 下 饱和 度 高 于 
通道 下 饱和 度 ， 扩 散 层 内 饱和 度 高 于 催化 层 。 同 样 ， 
有 接触 热 阻 时 ， 电 池内 温度 较 高 ， 所 以 ， 有 接触 热 
阻 时 ， 电 池内 液态 水 饱和 度 低 于 无 接触 热 阻 时 。 


现 


5 阴极 无 量 纲 饱 和 度 分 布 (a) 接触 热 阻 工 况 ，(b) 无 
接触 热 阻 工 况 


Fig.5 Water saturation distribution of cathode electrode (a) 


with TCR，(b) without TCR 
2.3 电池 极 化 曲线 
两 种 工 况 下 计算 得 到 的 电池 极 化 曲线 如 图 6 所 
示 。 电 压 高 于 0.6V 时 ， 两 种 工 况 极 化 曲线 无 差别 ， 
电压 高 于 0.6V 时 ， 有 接触 热 阻 时 电池 性 能 优 于 无 接 
触 热 阻 时 ， 当 电池 电压 为 0.3V 时 ， 二 者 相差 8%。 


| 浊 


一 一 有 接触 热 阴 
0 -0 - 无 接触 
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流 密 度 /Aecm” 

6 电池 极 化 曲线 (a) 接触 热 阻 工 况 ，(b) 无 接触 热 阻 工 况 
Fig.6 Cell polarization curve (a) with TCR , (b) without TCR 
值得 指出 的 是 ， 本 文 模型 未 涉及 接触 电阻 ， 因 
此 ， 极 化 曲线 的 差别 仅 由 温度 分 布 不 同 引 起 。 之 前 
的 分 析 指 出 ， 接 触 热 阻 使 电池 内 温度 升 高 ， 液 态 水 
饱和 度 降低 ,减轻 了 水 淹 现 象 ， 有 利于 反应 气体 的 
扩散 。 电 池 电 压 较 高 时 ， 浓 差 极 化 不 明显 ， 因 此 ， 
二 者 极 化 曲线 无 差别 ， 电 压 较 低 时 ， 浓 差 极 化 作用 
逐渐 加 强 ， 有 接触 热 阻 工 况 饱 和 度 低 ， 有 利于 传 质 ， 
因此 性 能 较 无 接触 热 阻 时 好 。 


现 


3 结论 

本 文采 用 三 维 两 相 非 等 温 燃 料 电池 计算 模型 分 
析 了 接触 热 阻 对 电池 性 能 的 影响 ， 计 算 结 果 表 明 : 
(1) 接触 热 阻 对 电池 内 温度 分 布 有 很 大 影响 ， 有 接 
触 热 阻 工 况 ， 电 池 最 高 温度 为 36SK， 电 池 肋 板 内 温 
度 几 乎 保持 不 变 ，MEA 内 温度 近似 为 ”型 分 布 ; 
而 无 接触 热 阻 工 况 ， 电 池 最 高 温度 为 336.SK， 电 池 
肋 板 内 温度 从 边界 开始 逐渐 升 高 .MEA 内 温度 星 “A” 
型 分 布 。 因 此 ， 模 拟 温 度 分 布 时 不 能 忽略 接触 热 阻 
影响。 
(2) 有 接触 热 阻 工 况 ， 电 极 液态 水 饱和 度 较 低 ， 此 
外 ， 由 于 模型 未 涉及 接触 电阻 ， 电 池 性 能 略 高 于 无 
接触 热 阻 工 况 ， 最 大 差别 约 8% 左 右 。 
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